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ABSTRACT 

We report on an analysis of the gas and dust budget in the the interstehar medium 
. (ISM) of the Large Magellanic Cloud (LMC). Recent observations from the Spitzer 

' Space Telescope enable us to study the mid-infrared dust excess of asymptotic giant 

0^ , branch (AGB) stars in the LMC. This is the first time we can quantitatively assess the 

■ gas and dust input from AGB stars over a complete galaxy, fully based on observations. 

The integrated mass-loss rate over all intermediate and high mass-loss rate carbon-rich 
AGB candidates in the LMC is 8.5 x IQ-^ MQyi-\ up to 2.1 x IQ-'^M^yr-^. This 
number could be increased up to 2.7 x 10~^ Mq yr~^ if oxygen-rich stars are included. 
5^ I This is overall consistent with theoretical expectations, considering the star formation 

rate when these low- and intermediate-mass stars where formed, and the initial mass 
functions. 

AGB stars are one of the most important gas sources in the LMC, with supernovae 
(SNe), which produces about 2-4x 10~^ yr~^. At the moment, the star formation 
rate exceeds the gas feedback from AGB stars and SNe in the LMC, and the current 
star formation depends on gas already present in the ISM. This suggests that as the 
gas in the ISM is exhausted, the star formation rate will eventually decline in the 
LMC, unless gas is supplied externally. 

Our estimates suggest 'a missing dust-mass problem' in the LMC, which is simi- 
larly found in high-z galaxies: the accumulated dust mass from AGB stars and possibly 
SNe over the dust life time (400-800 Myrs) is significant less than the dust mass in 
the ISM. Another dust source is required, possibly related to star-forming regions. 



Key virords: galaxies: evolution - galaxies: individual: the Magellanic Clouds - (ISM:) 
dust, extinction - stars: AGB and post-AGB - stars:mass-loss - (stars:) supernovae: 
general 
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1 INTRODUCTION 

The interstellar medium (ISM) of a galaxy is one of 
the drivers of its evolution, and the composition of the 
ISM determines many of the characteristics of the next 
generation of stars. The ISM is itself continuously re- 
newed and enriched by stellar ejecta. The enrichment 
occurs as stars evolve and die, either exploding as su- 
pernovae (SNe) or experiencing intense mass loss in a 
super-wind. Super-winds occur in low and intermediate 
mass stars during the asymptotic giant branch (AGB) 
phase (main sequence masses in the approximate range 
1-8 Mq), and in more massive red supergiants. The 
ejecta are enriched with elements produced in various 
phases of nuclear burning. In general terms, gas ejected 
from SNe includes newly synt hesised heavy elements , 
such as oxygen, iron and silicon (jNakamura et al.iri999f ). 
while AGB stars (below SM©) synthesize lighter ele- 
ments, especially carbon and nitrogen (Macdcr 1992). 
The chemical evolution of the gas can be well un- 
derstood in terms of the different stellar sources (e.g. 
IChiappini et al.l[200lh . taking into account the need for 
infa ll of unproces sed gas to stabilize the final metallic- 
ity ([Finlator fc Davee..2008) . The chemical evolution of 
the Large Magellanic Cloud (L MC) has recently been 
studied bv lCarrera et all (|2008( l. 

The origin of the dust in the ISM is less well un- 
derstood. Dust forms in stellar ejecta at temperatures 
of about 1000~1500K and at high densities. Different 
stars produce different types of dust. Red supergiants 
produce oxygen-rich (silicate) dust. AGB stars have two 
distinct chemical types: oxygen- and carbon-rich, de- 
pending on the abundance ratio of oxygen and carbon 
atoms within their atmospheres. Oxygen-rich stars form 
silicate dust, while carbon-rich stars yield amorphous 
carbon, graphite and SiC dust. Supernovae can produce 
both dust types, depending on the abundan ces in the 
different layers of the ejecta (|Rho et al.ll2008l ). 

The current rate of ISM enrichment by dust and 
gas depends on the total stellar population, the ini- 
tial mass func tion and the star-forming history (e.g. 
ISalpeteilll955| ). Type II SNe are expected to dominate 
the e nrichment in th e early phases of gal a xy ev olution 
(e.g. lMaedeij[l99^ . iHiras'hita & Ferrara| (|2002[ ) argue 
that dust grains injected from SNe into the ISM accel- 
erate star formation in young galaxies. It takes more 
than 100 Myr for the first intermediate-mass stars t o 
evolve onto the AGB fe.g. IVassiliadis fc Woodlll993D . 
Thus, dust and gas enrichment from AGB stars occurs 
later than from high-mass stars. Different galaxies, at 
different stages of this process, may be expected to show 
differences in gas-to-dust ratios, dust content, and, in 
consequence, ISM dust extinction curves. 

In our Galaxy, the major dust sources are presumed 
to be AGB stars and SNe (|Gehrzl [l989l ). Some other 
sources, such as Wolf-Rayet stars and novae, also con- 
tribute dust to the ISM of the Milky Way, but only 
in small quantities (|Gehr j [l989[ ). The relative i mpor- 
tance of AGB stars and SNe remains uncertain. iDwekl 
()1998( ) suggests that SNe are important silicate dust 
sources, while most car bonaceous dust gra i ns are from 
carbon- rich AGB stars. [Jura fc KleinmannI (jl989[ ) show 



that AGB stars are an important gas source within our 
Galaxy, except in the Galactic plane where SNe appear 
to dominate. The dust fo rmation efficiency in S Ne re- 
mains controversial (e.g. ISugerman et aL 2006t ). since 
SN shocks also destroy dust ( TielensI 1994 ) and since 
dust production prior the explosion remains unclear. 

There are also uncertainties in the dust formation 
efficiencies for AGB stars. This is expected to depend 
on the mass-loss rate, chemical type, and metallicity. 
For oxygen-rich st ars mass-loss rate s decrease towards 
lower metaUicity (IWood et al.l ll998t H owen fc WillsonI 
1991!; Marsh gdl et al.l 12004 ). In contrast, for carbon- 
rich stars, although the metallicity dep e ndence of mass- 
loss rates is still unclear (lHabinelll996l:lvan Loonl l2000l: 
IWachter et al] 120081: iMattsson et al.ll2008D . there does 
not appear to be any obvious metallicity dependence in 
the range fr om solar, down to one- twentieth of the solar 
meta l licitv (iGroe newegcn et al.l 120071 : iMatsuura et al.l 
l2007t ISloan et anr2009. ). However, we would expect the 
mass-loss rates of carbon stars to differ between envi- 
ronments with different abundances of alpha elements. 
In these cases the initial abundance of oxygen in the 
carbon-star precursor would be vital in determining 
how much carbon must be dredged up before any was 
free to become dust. 

Galaxy evolution models aim at fitting the current 
composition of the ISM by following the evolution of the 
stellar population over the lifetime of the galaxy. This 
requires a detailed knowledge of dust formation efficien- 
cies, gas and dust expulsion efficiencies and nucleosyn- 
thesis in stellar interiors. Here we aim to measure the 
current integrated rate of gas and dust input into the 
ISM for one specific galaxy: the LMC. Recent observa- 
tion s by the Spitzer Space Telescope (hereafter Spitzer 
(Wer ner et al.|[20M )) facilitate the detection of individ- 
ual dust formation sites in nearby galaxies. Observa- 
tions from Spitzer can thus be used to provide a new 
estimate of the gas and dust budgets for these galaxies. 

We have obtained spectra of AGB stars (mainly 
carbon-rich stars) in the LMC and the Small Mag- 
ellanic Cl£ud_(SMC}_using the Infrared S pectrograph 



ellani c t^ioud (blVlU) u smg tne intrared b pectrograpn 
(IRS ; iHouck et al.ll2004l) o n-boa rd Spitzer jsioan et ahl 
20061 120081: IZiilstra et all l2006t [Matsuura et al.l 12006 



Wood et al II2007I : iLagadec et al.H2007l: iLeisenring et al! 



2008[) : there are also independent spectral surveys 
of th ese galaxies (| Buchanan et al.l 120061 : iKastner et al.l 
I2008D . A complete LMC photo metric survey using 
Spitzer was recently published by iMeixner et al.l (j2006l 
SAGE). Combining these data we obtain a census of 
the mass-losing stars, and can then compare AGB stars 
with other sources of dust and gas in the LMC. The 
SMC data are not discussed in detail, but are used to 
extend the colour mass-loss relations to lower mass-loss 
rates. 



2 ANALYSIS 

2.1 Mass-loss rates and colours 

It has been shown empirically that mass-loss rates of 
AGB stars correlate with their infrared colours (e.g. 
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Figu re 1. Gas mass-loss rates for L MC and SMC c a rbon stars 
from ICroeneweeen et al. I ||2007| ') and ICruendl et al. I ||2008|') are 
lotted as a functioii of [3.6] — [8.0] colours extracted from SAGE 
Meixner et"alll2006l) . A gas-to-dust ratio of 200 is adopted. The 
solid (red) curve is the fit to the LMC sample only and the dashed 
(green) curve is the fit to the combined LMC and the SMC sam- 
ples. 
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Figure 2. Similar to Fig.[T] showing [3.6] - [24] for the LMC 
sample. 



l,Tiiralll987tlWhitelock et allll994l : lLeBertrelll997t ). Par- 
ticularly well known is the correlation with K — [12] 
(|Whitelock et all [20061 : | van Loon et al.1 I1999D . where 
[12] is the IRAS 12 /xm magnitude. It is assumed that 
this arises because K magnitude is a measure of the 
flux from the central star, while [12] measures the 
flux from the circumstellar envelope. A similar mass- 
loss rate versus colour relation has been obtained us- 
ing [6.4] — [9. 3] for carbon-rich AGB st ars in the LMC 
and S IVIC (jGroenewegen et all [2003). ICroenewegenl 
((200I) also provide theoretical relations between mass- 
loss rates and in frared colours for pho tometry from 
Spitzer. 2MASS (|Sknitskie et al.l l2006h and AKARI 
(jlVlurakami et al.ll2007[ ). 

We derive observational mass-loss rate versus colour 
relations by adopting gas mass-loss rates estimated 
for ca rbon-rich AGB sta rs in th e IVlag ellanic Clouds 
from iGroenewegen et al] (|2007D and iGruendl et al.l 
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Figure 3. Similar to Fig.[T] showing Ks — [8.0]. 



("20081). IGroenewegen et al.l (|2007t) used the spectral 
ener gy distribut ions obtained from Spitz er IRS spec- 
tra (jSloan et al.ll2006l : IZiilstra et al.ll2006D and JHKL 
mags observed at Siding Spring Observatory (SSO) in 
Austr alia, near simulta neously with the IRS observa- 
tions ([Wood et al.ll2007l ). The spectral energy distribu- 
tions were fitted using a radiative transfer model. The 
gas-to-dust ratio is assumed to be 200 for these carbon- 
rich A GB stars. For all these stars, we ext ract Spitzer 
IRAC (jFazio et al.ll200l and MIPS (Rieke et al.'!20Ql 
mags. For the LMC the Spitzer and the 2MASS mags 
are taken from the 2007 version of SAGE (Surveying th e 
Agents of a Galaxy's Evolution; iMeixner et al.ll2006D. 
while for the SMC they are from S^MC (|Bolatto et al l 
[2007'). Stars in clusters are excluded here so as to 
avoid possible source confusion. We cross-identified for 
sources within lOarcsec from SAGE and SfMC, us- 
ing the coord i nates quoted in ISloan et"al] (|2006D and 
IZiilstra et all ()2006t) . We found 24 stars in the LMC 
and 19 in the SMC. Some of them have incomplete pho- 
tometric data sets from 2MASS and/or SAGE. Highly 
reddened stars were n ot detected in the 2MASS J-band. 
iGruendl et al.l (|2008[ ) discovered thirteen extremely red 



objects (EROs) in the LMC, and seven of them were 
confirmed to be carbon- rich AGB stars, based on IRS 
spectra; we use their mass-los s rates in the following 
analysis. IGruendl et al.l (|2008f ) measured their magni- 
tudes in the Spitzer IRAC bands and the MIPS 24 and 
70 fim bands, but these stars we re not detected in the 
near- infrared either by 2MASS (jSkrutskie et al.ll2006l) 
or bv lKato et al] (|2007t) . 

The average interstellar extinction is estimated to 
be E{B - F)=0.15 and Aks^OM mag t owards the 
LMC and Aks=OM mag towards th e SMC ( Glass et al] 
ll999HZaritskvll999l : ICioni et al.l2006a.b ). We can safely 
ignore the interstellar extinction here, except for the 
J — Ks versus Ks diagram in Fig. lAll for a comparison 
with a figure in[Cioni et al. (2006a). 

Figs. [U - [3| show the mass-loss rates as a func- 
tion of [3.6] -[8.0], [3.6] -[24] and ifs- [8.0], respectively. 
These relations are fi t ted w ith a log function, following 
iLeBertre fc WintersI (|1998[ ). and the fits are shown in 
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the figures. Fits to the LMC and SMC combined sample 
give 

log Mg = -6.20/(([3.6] - [8.0]) + 0.83) - 3.39 (1) 
in the range of 1 < [3.6] - [8.0] < 9 



log Mg = -14.50/((i4:s - [8.0]) + 3.86) - 3.62 



(2) 



in the range of 1 < if s - [8.0] < 9, and fits to the LMC 
sample (3 < [3.6] - [24] < 14) give 

log Mg = -9.89/(([3.6] - [24]) + 1.43) - 3.36 (3) 

where Mg is the gas mass-loss rate. These equations are 
fitted for stars with colour ranges of 1.5 < [3.6] — [8.0] < 
9 for eq.C]) 2 < - [8.0] < 8 mag for eq.@, 2.5 < 
[3.6] - [24] < 13 for eq.®, respectively. All fits have 
correlation coefficients higher than 0.8. Coefficients are 
higher than 0.96 for [3.6] - [8.0] and [3.6] - [24] . The pho- 
tometric data of the EROs contribute to the better fit in 
the reddest colour range for the Spitzer bands, but these 
stars were not detected in Ks. MIPS [24] data are used 
for the LMC only. Relatively low mass-loss rate stars are 
not detected in MIPS [24] (Fig.[X5]), and the equation 
using [3.6] — [24] is suitable for relatively high mass-loss 
rate stars only. Stars in the SMC sample have bluer 
infrared colours and lower mass-loss rates than those 
of the LMC sample; this bias was present in the origi- 
nal sample selection for IRS spectroscopy (ISloan et aTl 
120061: IZiilstra et al.ll2006t iLagadec et all 120071) . If data 
in all bands are available the equation in [3.6] — [8.0] is 
the better one for deriving mass-loss rates, except for 
low mass-loss rate stars, where Ks — [8.0] will be more 
sensitive to the small dust excess. These two equations 
will be used later. 

These fits should be regarded as very uncertain 
in the extrapolated regime of extremely low mass-loss 
rates. Also, the equation fitted for Ks — [8.0] is not 
appropriate for the mass-loss rates of red carbon stars 
(> lO-^Moyr-i). 

Throughout this paper we assume a gas-to-dust ra- 
tio of 200 for carbon-rich AGB stars to convert the mea- 
sured dust mass to a gas mass-loss ra te. This ratio is 
simil a r to Galactic values (c.f. Habin g |1996t iLe Bertrd 
119971: ISkinner et £d]|l999[ ). [Matsuura et all (|2005D ar- 
gue that, for carbon-rich AGB stars, the gas-to-dust 
ratio only marginally depends on the metallicity of 
the host galaxy. This is because carbon atoms are 
synthesised in AGB stars, and the quantity of amor- 
phous carbon (or graphite) produced depends on the 
number of excess carbon atoms in the atmosphere of 
the AGB star after all the oxygen atoms are locked 
into the stable carbon monoxide molecule. Various ob- 
servations suggest no obvious differences in mass-loss 
rates for carbon-rich AGB stars amongst the LMC, 
the SMC and the Fornax dwarf Spheroidal G alaxy 
(jGroenewegen et all l2007t iMatsuura et all I2007D : the 
time variations and uncertainties of the mass-loss rates 
appear to be larger than any metallicity effects. Typ- 
ical metallicities of the LMC, the SMC, and the For- 
nax dwarf Spheroidal Galaxy are half, one quarter and 
one tenth of t he solar metallic i ty, res pect ively. Theoreti- 
cal models bv lMattsson et al.] (|2008l ) and lWachter et all 



||2008( ) predict that the mass-loss rates should depend 
on metallicity, although the predicted difference be- 
tween the LMC and SMC amoun ts to on l y a f actor 
of two. iHabind ([19961 and M atsuura et all |2007t) and 
argued that the mass-loss rate of carbon stars is pri- 
marily determined by the difference in the numbers of 
carbon and oxygen atoms: as the excess carbon is self 
produced by the AGB star, carbon stars will show rela- 
tively little difference in mass-loss rates between dif- 
ferent galaxies, although for low metallicity systems 
the mas s loss may begin earlier on the AGB phase 
(Lagadcc fc Ziilstral l2008al ). Note that this will apply 
only to carbon-rich AGB stars; the mass loss of oxygen- 
rich AGB st ars is like ly to sh ow a strong metaUicity 
dependence (IWood et a l. 1998: IBowen fc Willson|[ 19911 : 
Ivan Loon|[2"000l : ILagadec fc Ziilstrall2008a[ ). 



In the process of estimating the mass- loss rate 
(jGroenewegen et al.|[2007[: iGruendl et al.ll2008f) . an out- 
flow velocity of the circumstellar envelope is assumed. 
For oxygen-r i ch stars, lower outflow veloc ities are found 
([Wood et al.l [19921 : [Marshall et"al[ I2004D . but the gas 
velocity has never measured for carb on-rich stars in 
the LMC. Theoretic al models (Mattsson et al.[ [2008[ : 
[Wachter et al.l [2008f) suggest that the mean veloc- 
ity d oes not depend on meta llicity for carbon-rich 
stars. [Groenewegen et ap ([2007[) assumed a g as veloc- 
ity of lOkms"^ while [Gruendl et all ([20081 ) used 8- 
12kms"^ 

The other uncertainty originates from the as- 
sumption that the dust con d ensation temperatur e 
is 1000 K ((G roenew eee^ 120061: IGruendl et al] I2008D . 
This is a siinilar v alue as used for Galactic stars 
(Groene wegenI [1998I ). However, for example, the con- 
densation temperature of graphite dust grains varies 
with carbon-to-oxygen (C/0) ratio ([Lodders fc FeglevI 
[19951: [Tielens et al] [2005[ ). and at a C/0 ratio of 1.2, 
the temperature is 1700 K. It has been suggested that 
the C/0 ratio of carbon-rich stars in the LMC is sys- 
tematically higher than in t heir Galactic counterparts 
([Matsuura et al J [20051 [200I) . This has also been sug- 
gested fro m the observations of PNe in the LMC and 
the SMC (iLeisv fc Denn 'ifeIdl[T99l i Dooita et al.l[l997[ : 
|Stanghellim_eLalJl2005'). To account for C/0 ratio vari- 
ations and to test their influence of the dust condensa- 
tion temperature on mass-loss rate, we tested radiative 
transfer modelling of the SED fitting, with dust con- 
densation temperatured of 1000 K, 1200 K, 1300 K and 
1700 K. We used the spectrum of one of the highest 
mass-loss rate stars, 05(3231.49 -680535.8 (J2000), ob- 
served by [Gruendl et al.1 ([2008D . The result shows that 
both the fitted optical depth and expansion velocity in- 
crease with dust condensation temperature. The mass- 
loss rate required to fit the SED could be a factor of 2.4 
higher, if the dust condensation temperature increases 
from 1000 K to 1700 K, although the numerical error 
is large for the models with dust condensation temper- 
ature above 1300 K. This suggests that the mass-loss 
rate, which is estimated from infrared colour, could be 
underestimated by a factor of up to 2.4, systematically. 
In the overall analysis, we give the results under the 
assumption of 1000 K dust condensation temperature. 
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however, in the final gas and dust budget, we include 
this potential increase. 

AGB stars are long period variables, whose magni- 
tudes vary over a 100-1000 day time scale. The ampli- 
tudes of the infrared variability are smaller than those 
of optical variabilitie s: they range from 0.2 to 1.5 mag 
at L -band (3 fim) ()Le Bertrd Il992t IWhitelock et all 
12006 ') . Furthermore, mass- loss ra tes of AGB stars vary 
over time scale s of 1000 years (jlzumiura et al.l 119961 : 
iLeao et al.|[2006l) . The cause o f long-term mass-los s rate 
variations are still unknown (jZiilstra et al.|[2003 ). The 
variability affects our analysis for individual sourcs, but 
are averaged in the full, large sample. 



126 red stars in the LMC, and identified 11 carbon- 
rich and 61 oxygen-rich stars. These are described as 
'C-AGB (optical)' and '0-AGB (optical)' in the fig- 
ures. We also included M-type stars (both giants and 
superg i ants) from fol l owing works, [S andulcak & Phili 
(119771). iBlanco et al.l (fl978h. IWesterlund et al.. a98lf 



IWood et al.l (|198l I1985D . iHughed (I1989D . iReid et alJ 
(| 19881 ) using coordinates taken from SIMBAD. There 
are 1710 objects, having 1568 counterparts in IRAC and 
971 counterparts in MIPS within 1 arcsec. These stars 
are identified as '0-AGB/RGB (Sim)' in the figures, al- 
though some of them could be red giant branch (RGB) , 
rather than AGB stars. 



2.2 Separation of oxygen-rich and carbon-rich 
AGB stars 

To estimate the integrated mass-loss rates from carbon- 
rich AGB stars, it is essential to find a classification 
scheme that will separate carbon-rich from oxygen-rich 
stars. In particular, oxygen-rich AGB stars and red 
supergiants (RSGs) follow different r aass-loss rate vs. 
colou r relations from carbon-rich stars (jWhitelock et al.l 
|2006() and their dust contents will be very different. 

We first compiled a list of AGB stars in the 
LMC observed with the Spitzer IRS by various groups 
(IZiilstra et all 120061: ISpeck et al.l [20061 iKastner et all 



120081: iLeisenring et al.ri2008l - ISloan et al.l l2008l Kem- 
per et al. in preparation). These studies selected stars 
mainly from previo us LMC mid-infr ared surveys, IRA S 
(ISchwering fc Israeli Il990t ). MSX (lEgan et all l2001l) 
and s urveys for long-period variables fe.g. lHugheslll989l 
Il990t ). We used the carbon- and oxygen-rich classifica- 
tions quoted in these works. These are mostly based 
on the dust features (either silicate or SiC at 10 /xm), 
but for stars without dust features, molecular bands are 
used for chemical classifications. The group of oxygen- 
rich stars includes both giants and super-giants, but 
as shown in Sect lAl) two young stellar objects (YSOs) 
or post-AGB stars with silicate bands may contami- 
nate the sample. We searched the photometric data in 
SAGE to find the counterparts to the stars observed 
with IRS: 122 stars have counterparts within 1 arcsec, 
7 stars within 1-2 arcsec, and 5 stars have more distant 
counterparts. We use only identifications where the po- 
sitions agree to within 2 arcsec. For the LMC this pro- 
vides a total of 76 oxygen-rich and 40 carbon-rich stars 
from the SAGE database. We further add one oxygen- 
rich star (IR AS 05298—6957) fro m an OH maser survey 
in the LMC (jWood et al.lll992D : the other known OH 
maser sources were observed with IRS. In the following 
figures the OH/IR star is included in '0-AGB/RGB 
(IRS)'. 

A catalogue of spect roscopically identifie d optical 
carbon stars in the LMC ([Kontizas et al.l[200l[ ) was also 
searched for SAGE sources within 3 arcsec. In practice, 
98 percent of the sources are found within 1 arcsec, 
yielding 5710 SAGE counterparts. Where there are 
more than two SAG E sources within 3 arcsec of the 
iKontizas et al.l (|200lf ) coordinates, we chose the closest 
one.^ 

ICioni et al.l (j200ll ) have obtained optical spectra of 



2.2.1 Classifications involving the 8 band 

In this section we develop a classification scheme us- 
ing the [3.6] — [8.0] colour to separate carbon-rich from 
oxygen-rich stars. Diagrams involving other colours can 
be found in the Appendix. Fig. 0] shows a [3.6] — [8.0] vs 
[8.0] colour magnitude diagram. The [8.0]-magnitudes 
of carbon-rich stars peak at [3.6] — [8.0] =2-3 mag. The 
SED of thick circumstellar shells peaks longward of 
about 7 /xm. 

In FigH] oxygen-rich stars tend to show brighter 
[8.0]-mags than carbon-rich stars at any given 
[3.6] — [8.0], with the exceptio n of the two oxygen - 
rich post-AGB/YS O candidates. [Meixner et all (|2006D : 
iBlum et all (|2006f) used the [3.6] - [8.0] vs [8.0] dia- 
gram for sou rce classification. W e also classify 'extreme- 
AGB stars' ([bI um "eTall [20061) as oxygen- or carbon- 
rich, based on our Spitzer IRS spectra. The boundaries 
which we use to separate carbon- from oxygen-rich AGB 
stars in Figd] are: 

[8.0] > 9.0 - 1.8 X [3.6] for 0.7 < [3.6] < 2.5 (4) 
[8.0] < 11.2- 1.8 X [3.6] for 0.7 < [3.6] < 2.0 (5) 
[8.0] > 4.0 - 0.2 X [3.6] for 2.0 > [3.6] (6) 
and 

[8.0] < 7.2 - 0.2 X [3.6] for 2.5 > [3.6]. (7) 

This region is in general consis tent with the evolu tion- 
ary tracks of carbon-rich stars (Ma rigo et al.|[2008l ). 

In addition to low luminosity oxygen-rich stars 
fSect. ETj) . nine IRS classified oxygen-rich stars are 
found in the carbon- rich re gion. Four of them are known 
to be OH/IR stars (Woo d et al]|1992D and their IRS 
spectra show the silicate band in absorption. It is not 
possible to make a perfect division of the sources by this 
method, and inevitably some oxygen-rich sources are 
found in the carbon-rich region. However, this colour 
selection provides a method of isolating the majority of 
carbon-rich stars. We note that the high mass-loss rate 
OH/IR stars found in the carbon rich region are not ex- 
pected to occur in large numbers. Furthermore, carbon- 
rich stars with colours similar to the OH/IR stars have 
lower mass-loss rates than the OH/IR stars, due to their 
different dust properties. Thus, the influence of blended 
oxygen-rich stars would be only a small perturbation 
on our estimate of the integrated mass-loss rate from 
carbon stars. 
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Figure 4. [3.6] — [8.0] vs [8.0] colour-magnitude diagram. The top panel shows the colour magnitude diagram for the complete SAGE 
sample. In the lower panel, spectroscopically identified oxygen-rich and carbon-rich AGB stars are plo tted. The selection cr iteria used to 
identify carbon-rich AGB candidates are indicated by the solid lines. The sensitivity limits quoted bv lMeixner et a~ for the final 

data product are indicated by the dashed lines. [8.0] ~ 10.5 likely corresponds to the tip of the RGB branch the number of carbon stars 
decreases at about [3.6] — [8.0]~0.7 (at about [8.0]~9), As seen in the top panel, 
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There are a large number of so called 'optical' 
carbon-rich stars outside the region defined by equa- 
tions (jUIll). However, with the available colour-colour 
or colour-magnitude diagrams we are unable to sepa- 
rate oxygen- and carbon-rich stars with thin circum- 
stellar envelopes. Mass loss from carbon-rich stars with 
thin shells is discussed separately fSect. [5TT|) . 

At around [3.6] — [8.0] > 4, a small number of un- 
resol ved galaxies might intrude into the AGB star re- 
gion. iBlum et all (|2006D investigated the [3.6] - [8.0] 
colours of distant galaxies. They found galaxy candi- 
dates in the red and low luminosity regime, concentrat- 
ing at [3.6] - [8.0] ~ 3 and [8.0] > 10 (their fig. 7). 
They suggested that galaxy candidates extend up to 
[8.0] ~ 7 at [3.6] - [8.0] - 3. They did not classify 
point sources with [3.6] — [8.0] ^ 3 and [8.0] < 10. How- 
ever, our sample based on Spitzer IRS spectra clearly 
shows the presence of carbon-rich stars at [8.0] ~ 7 at 
[3.6] — [8.0] ~ 3. We argue in Sec l3.2l that the majority 
of objects with [8.0] < 7 and [3.6] - [8.0] > 3 belong 
to the LMC. From this we might also deduce that the 
most populous bright near- and mid-infrared sources in 
galaxies are AGB stars. This will im pact the infrared 
colour of galaxies (jTonini et al.ll2008l ). 

There is a change in the number density of car- 
bon stars at Ks— [8.0]~1.4 (top panel of Fig. lA3[) and 
[3.6] - [8.0] ~ 0.7 (top panel of Fig.©. This blue/red 
colour separation might correspond to AGB stars with 
and without significant mass-loss rates. For the bluer 
sources the opacity includes molecular bands and H^ 
continuous absorption. For the redder sources, dust 
emission contributes to the mid-infrared excess. In- 
deed, IRAS observatioi is show a [12] — [25] gap be- 
tween -1.2 and -1.4 (jvan der Veen fc Habind I1988D 
and stars redder than —1.4 are those with dust-driven 
mass losffl The gap can be understood as the sepa- 
ration of huge numbers of stars with very little mass 
loss from those with circumstellar envelopes. Once dust 
grains have condensed, m ass-loss rates increase drasti- 
cally (Ivan der Veen fc Ha bing 1988), as do colours. The 
Ks — [8.0] and [3.6] — [8.0] do not show the obvious gap 
seen in [12] — [25], as the influence of dust on the [8.0]- 
band is not strong, but it is likely that the decrease 
in stellar density is related to the increase in mass-loss 
rate. 

Possible sources of contamination f or our carbon- 
rich A GB sample include Hll regions. iKastner et al.l 
()2008D have classified bright objects in the MSX A- 
band, which covers 8 /xm, and they found four Hll re- 
gions in the range 4.0 < A < 6.5 and 6 < K — A < 9. 
Although the MSX A- and Spitzer [8.0] -bands are not 
identical, the wavelength coverage is similar, suggesting 
Hll regions could be a contaminant. A search for cross- 
identifications of SIMBAD sources within three arcsec 
revealed one Hll region, five Be stars or emission-line 
objects, two YSOs, three molecular clouds, one sym- 



^ In the van der Veen et al. study the IRAS bands are not 
normalized with zero-magnitudes, but [12] — [25] is defined by 
2.51og(F25/-Fi2), where F12 and F2S are the IRAS fluxes in Jy in 
the 12 and 25 ^m bands. 



biotic star (essentially an AGB star with a compact 
com panion), two planetary nebulae (PNe), and one WR 
star. lWhitnev et all ()2008f ) showed that luminous YSOs 
cross the AGB sequence, but we expect the high mass 
YSO population to be significantly smaller than that of 
AGB stars. The initial masses of carbon- rich AG B stars 
are about 1.5-5.0Mq (jVassiliadis fc Woodl[T993l) . while 
YSOs, which can reach [8.0] < 8, have m asses higher 
than 10 Mq. Using the IMF from .Kroupal (|200L ). with 
an exponent of a = 1.3, and assuming a constant star 
formation rate, stars with 1.5-5.0 have a birth rate 
ten times higher than stars with ma sses greater than 
10 M©, although the SFR might vary (|Feastlll995[ ). The 
lifetime of such high mass YSO is also much shorter 
than that of the AGB phase. Similarly WR stars have 
high initial-masses and their numbers should be low. 
PNe (and most symbiotic stars) are usually located 
outside of our carbon-rich region, and are more con- 
centrated at lower [8.0]-band luminosities (|Hora et al.l 
200^ who used the m ost complete PN catalogue from 
[Reid fc Parkeil (|2006f )). We expect that because their 
lifetimes are short, the PN contamination is low. 

The sensitivity limits of the surveys are shown on 
the colour-magnitude diagrams. These sensitivities (5tT 
detec tion limits) are those for t he final produc t take n 
from lSkrutskie et all (L2006) and 'M eixner et all (|2006D : 
note that the actual detection limits of the catalogues 
used could be lower than those illustrated here. The 
Spitzer SAGE project includes observations of the LMC 
at two epochs, but the catalogue used (2007 Winter 
vers ion) contains on ly data from the first of epoch 
(■Meixner et al.ll2006D . In Fig.EH the sensitivity (5g 
detec tion) limits of the 2MASS survey (jSkrutskie et al.l 
120061) are plotted on the diagram. The 2MASS sensitiv- 
ity is not high enough to detect the reddest part of the 
carbon-rich sequence. 



3 RESULTS 

3.1 Mass-loss rates for all carbon-rich AGB 
candidates 

The mass-loss rates for carbon-rich stars are estimated 
using the colour vs mass-loss-rate relations derived in 
Sect. 12. II The stars with high mass-loss rates are listed 
in Table[T] We assume that statistically these relations 
will give an integrated mass-loss rate for all carbon- 
rich mass-losing AGB candidates. Specifically, we use 
[3.6] — [8.0], supplemented by Ks — [8.0] for low mass- 
loss rate stars, to estimate the rates. 

There are 1779 carbon-rich AGB candidates with 
[3.6] — [8.0] > 0.7 mag. The integrated mass-loss rate 
from all 1779 stars, derived from this colour, is 7.4 x 
lO~^M0yr~^. The highest individual mass-loss rate is 
9.3 X 10~^ Mq yr~^. The upper limit may be set by the 
detection limit in the [3.8]-band (Figs. |4| and E3)- The 
distribution of stars amongst the mass-loss rate ranges 
is summarized in TablejH 

For stars with low mass-loss rates, the fits to 
colour vs mass-loss rate using the two different colours, 
Ks - [8.0] and [3.6] - [8.0] > 0.7, give slightly differ- 
ent results; this is because towards blue colours, the 
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Table 1. LMC carbon-rich AGB candidates with high mass-loss rate (> 1 X 10~^ Mq yr~^). The full table is available in on-line. 
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Table 4. Overall properties of the LMC 



All surveyed area Bar Bar cent Disk 



Total number of carbon-rich AGB candidates 
Integrated mass- loss rate* ( X 10"'' Mq yr^^) 
Area (deg^) 

Mass-loss rate per unit area (xlO~* Mq yr~^ deg~^) 
Same as above but in different units (xlO"** Mq yr"'^ kpc~ 
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1.1 



The mass-loss rates from carbon-rich AGB stars could be underestimated by a factor of up to 2.4, due to the unknown dust 
condensation temperature (Sect l2.l"1l .t gas mass-loss rates are estimated using Ks — [8.0] and [3.6] — [8.0]. 



Table 2. Number distribution of mass-losing carbon-rich can- 
didates in the LMC. Two different colours {Ks — [8.0] and 
[3.6] — [8.0]) wore used to estimate the mass-loss rates of indi- 
vidual stars. In the last column, one or other of these colours 
were used, i.e., if [3.6] — [8.0] < 2.2, Ks — [8.0] vs mass-loss rate 
relation was adopted, otherwise the [3.6] — [8.0] vs mass-loss rate 
relation was used. The total gas mass-loss rate is obtained by 
integrating the mass-loss rates of individual stars. 
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Total 


1779 
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86§ 



A/gi gas mass-loss rates, assuming a gas-to-dust mass ratio of 200 
A^: number of stars 

A^gi : Mg estimated from [3.6] - [8.0] only 

Mg2 ■■ Mg estimated from Ks - [8.0] and [3.6] - [8.0] 

t The minimum mass-loss rate is 4 X 10~*AfQ yr~^, limited by 

[3.6] - [8.0] > 0.7 mag 

t The minimum mass-loss rate is 4 X 10~*AfQ yr~^, limited by 
Ks - [8.0] > 1.4mag (as well as [3.6] - [8.0] > 0.7mag). 
§ Due to rounding, adding Mg in seven ranges does not give the 
identical number as the total Mg. 



Table 3. Number distribution of LMC carbon-rich candidates as 
a function of mass-loss rates based on [3.6] — [8.0] and Ks — [8.0] 
colours in the entire SAGE survey area, the bar region, the bar 
centre, and the remainder of the surveyed area except of the bar 
region (quoted as 'disk') 
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-4 


0.7% 
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0.8% 


0.7% 


Total N 




1779 


1194 


517 


585 



Afg2 : gas mass-loss rates estimated using Ks — [8.0] and [3.6] — 
[8.0] 

N: number of stars 

Percentages over the seven ranges do not sum to 100 % due to 
rounding. 



extrapolated relations are poorly defined. Ks ~ [8.0] is 
probably better for these low mass-loss rate stars. We 
combine the two estimates of the mass-loss rate as fol- 
lows: We adopt the estimate from [3.6] — [8.0] for stars 
with [3.6] — [8.0] > 2.2, otherwise we take the rates 
from Ks — [8.0]. If a star does not have a Ks-ha,nd 
magnitude, we use the rates from [3.6] — [8.0]. Above 
[3.6] - [8.0] > 2.2 (- 3 X IQ-^M^yi--^) up to Ks de- 
tection limit, the mass- loss rates from the two colours 
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are essentially identical. This procedure gives an inte- 
grated mass-loss rate of 8.6 x lO~^M0yr~^, as sum- 
marized in the last column of Table[2l This combined 
method only makes a 10 % difference in the integrated 
mass-loss rate from that using the [3.6] — [8.0] relation 
alone, as the low mass-loss rate stars contribute little 
to the integrated value. 

Independent wo rk with a totally different approach 
bv lSrinivasan et al] ((2008) obtained a similar integrated 
mass-loss rate for AGB stars (2.5 x 10^^ M© yr""'^). Our 
estimate is higher, because we included stars with ex- 
tremely high m ass-loss rates, wh ich had not been cate- 
gorized before ( Blum et al.ll2006D, but were found to be 
carbon- rich AGB stars recently (|Gruendl et al.|[2008r ). 

Because of the sensitivity limit in the if s-band 
(Fig.EH, and possibly in the [3.6]-band (Fig.Ell, it 
is possible that not all AGB stars with extremely high 
mass-loss rates will have been detected. Therefore the 
integrated mass-loss rate derived above might be con- 
sidered a lower limit. 

Table[2] shows the number distribution per mass- 
loss rate range among carbon-rich AGB candidates, 
combining the two infrared colour estimates. Stars with 
relatively high mass-loss rates (> 6.0 x 10~® A/0yr~^) 
are responsible for most of the integrated mass 
loss. A similar result was foun d for Galactic AGB 
stars (iKnaPD fc Morrisl[r985: Jura fc Kleinmannlll989t 
iLeBertre et al.ll200ll: lOiha et al.ll2007i r 

There are 1779 carbon-rich AGB candidates in the 
LMC (SecEJ) selected by [3.6] - [8.0] > 0.7. The num- 
ber of AGB stars with [3.6] - [8.0] < 0.7 is almost 10 
times greater than that with [3.6] — [8.0] > 0.7. In the 
region of 10.5 < [8.0] < 8 and [3.6] - [8.0] < 0.7, there 
are 13460 stars, where we assume that [8.0] ^ 10.5 cor- 
responds to the tip of the RGB branch (as the number 
of stars decreases at [8.0] < 10.5; Fi^.j^, and where 
both oxygen-rich and carbon-rich stars are mixed to- 
gether. Furthermore, there are carbon-r ich AGB stars 
below the tip of the RGB ([8.0] < 10.5) (iLaga dec et aP 
I2007D . which are uncounted due to the confusion with 
RGB stars. These will be a mixture of extrinsic carbon- 
stars, which have become carbon-rich due to mass- 
transfer from an evolved companion, or AGB stars in 
the brief phase of low luminosity which follows a he- 
lium shell flash. Extrinsic carbon stars should not be 
included in our census. The numbers of the others 
could be significant, but their mass- loss rates are low. 
The boundary at [3.6] — [8.0] ^ 0.5 may be separat- 
ing stars with dust-driven winds from those with pul- 
sation/opacity driven mass loss. If we assume a plau- 
sible mass-loss rate of 1 x lO^^Moyr"^ on average, 
similar to red-giant branch stars (10~^-10~* Mq yr~^ ), 
the integrated mass-loss rate from these stars would be 
~ 1 X lO~'*M0yr~^, a minor contribution to the inte- 
grated mass-loss rate. These are very rough estimates, 
and we do not include them into the final number. 



3.2 Spatial distribution of carbon-rich stars 

Fig. [5] and Table[3] show the distribution of carbon-rich 
AGB stars with intermediate and high mass-loss rates 
within various parts of the LMC. The areas are specified 



6xl05/<Mg< 1x10* 




3xl(fi -f Mg < 6x10*-- 
" Mg < 1x10' 




Figure 5. The spatial distribution of mass-losing carbon-rich 
AGB candidates in the LMC is plotted over the SAGE 3.6 fim 
image (top and middle). Candidates are classified into groups ac- 
cording to mass-loss rate as given on the figure. Carbon-rich AGB 
candidates are concentrated in the LMC bar. In the bottom panel, 
large and small rectangles show the areas defined for the 'bar' and 
the 'bar centre', as used in our analysis, (on-line v ersion, the bot- 
tom p anel is the SAGE 3-color composition image llMeixner et al.l 
l2006h : 3.6 /im image for blue, 4.5 /^m image for green and 5.8 ^m 
for red. The LMC bar region is found in blue color, as stellar 
components are dominant.) 
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in the lower panel of Fig.O The majority of the stars 
are located towards the LMC bar. 

One-third of the mass-losing carbon-rich AGB can- 
didates are concentrated in the bar centre. This is con- 
sistent with the findings from near- infrared surveys, 
which sho w a high concentration of carbon-rich stars 
in the bar (|Blanco et al.lll978l: [Cioni et al.ll2006aD . 

There are at least two regions containing clusters, 
NGC2154 (RA=05h57m38s, Dec=-67°15.7', J2000), 
where about 44 ca rbon-rich stars have been reported 
(|Blanco et al.lll978D . and NGC 1978 (RA=05h27ml7s, 
Dec=-70°44.1', J2000 ) where about 24 carbon stars 
have been found by I Wester lund et al.l (|1991h . and 
which contains a carbon-rich AGB Mira variable 
(|Tanabe et al.lll997[ ). 



4 DISCUSSION 

4.1 Mass-loss rates from AGB stars per unit area 

Using the result listed in Table[4] for the LMC the inte- 
grated mass-loss rate is 2.2 x 10""* yr~^ kpc~^ for 
carbon-rich AGB candidates only. There is a spatial 
variation in the range 1.1-5.1 x lO""* yr~^ kpc^^. 
This value is compar able to that found in th e solar 
neighbourhood, where [jura fc KleinmannI (|l989t ) inves- 
tigated AGB stars within '^Ikpc from the Sun and 
obtained an integrated mass-loss rate from AGB stars 
of 3-6x10"'' Mgyr-^kpc-^. Their study was limited 
to stars with Mg> 2 x 10~®Moyr~^, but lower mass- 
loss rate stars would not contribute to the value signif- 
icantly. Half of the mass was from oxygen- rich and the 
other h alf from c arbon -rich AGB stars. .Thronson et al] 
(|1987D . iTielend ()l990f ) and iGroenewegen et all (|l99a 
obtained similar numbers for the gas injection rate 
from carbon-rich AGB stars in the Galaxy. Assum- 
ing that the gas-to-dust ratios are comparable in these 
two galaxies, the gas mass per unit area injected from 
carbon-rich stars in the Galactic plane and the LMC 
bar centre are similar. 



4.2 Other gas and dust sources in the LMC 

For comparison with the result discussed above, we esti- 
mate the gas and dust mass expelled from other sources 
in the LMC and summarize the findings in TablelHl 



4-2.1 Gas and dust injection from SNe 

We estimate the g as output from SNe based o n the the- 
oretical model of 'Kodama fc Ari niotol ()1997[ ). We use 
the revised IMF from Kroupa (200ll) (lower mass limit: 
Mio = 0.1 M0 and upper mass limit: M^p = 60 Mq), 
and assume a constant star formation rate (SFR) over 
the past few Myr which is the same as the current 
value. Under these conditions, the ratio of the SFR to 
the gas expulsion rate (G ER) by type II SNe is 0.32. 
ISmecker-Hane et al.l ()2002l ) measured the SFR using red 
giants, which date from as recently as 200 Myr ago; 
we adopt this as the current value. They gave SFRs 
of 0.01 Mq yr-i deg-^ and 0.004 M© yr"! deg-^ in the 



bar and disk regions, respectively. These numbers result 
in SN gas inputs to the ISM of 3 x 10"^ Mq yr^^ deg"^ 
in the bar and 1 x lO""* M© yr~^ deg~^ in the disk. This 
implies that, in the bar, the gas injection rate from 
SNe is higher than that from AGB stars, but AGB 
stars are more important in the disk, if oxygen-rich 
stars are counted. The fact that supernova remnants 
(SNRs) are concentrated ar ound the bar region and 
30 Dor (jWilliams et al.llT999f ) supports this conclusion. 

In the following we use the SN rate estimated from 
the number of SNRs and the dust mass observed from 
SNe and SNRs. Math ewson et al.l (1983) hst 25 SNRs 
in the LMC. They estimated the frequency of type II 
SNe as one per 275 yrs, and for type I SNe as approxi- 
mately one per 550 yrs. More recently, IWilliams et al.l 
(|1999( ) added a further 6 SNRs. Although their SN 
type is unknown, they are more likely to be Type 
II, because galaxy chemical evolution models suggest 
that Type la an d Type II occur with a ratio of about 
0.3 in the LMC (|Tsuiimoto et al.lll995D . iFilipovic et"al] 
(|1998l ) estimated an LMC SNR rate of one m every 
10 yrs. The Type II SN r ate (for both type I and II) 
of iMathewson et all (|1983D remains valid, but the time 
interval could be shorter by up to 145 yrs. 

The dust output from SNe remains uncertain. The 
famous type II SN in the LMC, SN 1987A, emit- 
ted a significant infrared excess from d ust grains 
(iMoselev et al.lll989l: IWhitelock et al.lll989[ ). The high- 
est estima te of the dust m ass from SN 1987A is 1.3 x 
IO-^Mq (jErcolano et all [2007). Based on this value, 
type II SNe could produce up to 5 x lO~^M0yr~^ 
of new dust. In the SMC dust mass of 8 x 1 0"^ M© 
is found in a SNR (|Stanimirovic et al.l I2005D . Even 
higher dust masses have been reported from other type 
II SNe in other galaxies: 0.02 M© in SN 2003gd in 
NGC 6 28 (ISugerman et af] 120061) and 0.02-0.05 Mq in 
Gas A mho et al.ll2008D . A lower dust mass has been 
sugge s ted for SN2003gd ( 4 x 10~^ Mq; iMeikle et al] 
HoOa)- ISakon et all (|2008l ) estimated a dust mass of 
2. 2-3.4X 10-3 Mq produced from SN 2006jc. These es- 
timates give a range of (0.1-130) xIO-^'Mq yr"^ for the 
dust mass from observations of SNe and SNRs. The dust 
mass produced by SNe remains very uncertain for the 
following reasons: extinction measurements give a col- 
umn density rather than a mass; the mid-infrared pho- 
tometric measurements are mainly sensitive to warm 
dust at the early phase, and may need to be cor- 
rected for flash- heated interstellar dust; sub- mm mea- 
surements probe cold dust, but can be confused with 
interstellar dust clouds. The clumping is also uncertain 
and affects the mass determination. Finally, the primary 
dust production by SNe may be significantly counter- 
acted by subs equent dust destruction when the ejecta 
hits the ISM (| Jones et al.l 11994 lDwekll2008l ). which is 
not considered here. The dust production given here 
for SNe is tentative, and could change by an order of 
magnitude or more. 

We adopt a lower li mit for the p rogenitor mass 
of Type II SNe_of 8 Mq (jSmMtLeLalj r2008) , and the 
IMF from iKroupal (|200l . with a GER/SFR through 
SNe of 0.32. Then, with a SN rate of one per 145-275 
yrs, the gas ejection rate from Type II SNe is 0.02- 
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Table 5. Gas and dust injected into the ISM of the LMC 



Sources 



Gas mass Dust mass 

(lO-3M0yr-i) (lO-6M0yr-l) 



Chemical type 



AGB stars 
Carbon-rich 
Oxygen-rich 

Type II SNe 

WR stars 

Red supergiants 



8.6 (up to 20.6§) 
»2t 
20-40 
0.6 

>1 



43 (up to 100) 
»0.4t 
0.1-130t both O- 

>2 



C-rich 
O-rich 
and C-rich 
(C-rich?) 
O-rich 



Considering the unknown dust condensation temperature, mass-loss rate can be underestimated by a factor of up to 2.4 (Sect l2.1 
t The number hstcd is the lower limit obtained from the observations. Expected gas and dust mass-loss rate are 
8.6 X 10-3 Mq yr-^and 1.2 X lO"'' Mq yr"!, respectively 
I Dust production (or possibly destruction) in and around SNe remains uncertain. 



O.O4M0yr~^. This is higher than the integrated gas 
mass-loss rate from carbon-rich AGB candidates. This 
number criticaUy depen ds on the p rogenitor mass range 
of SNe (see discussion in ISmartt et aL 2008). Type Ib/c 
SNe may also contribute to core-collapse SNe and are 
not counted here as their origin remains uncertain. 



4.2.2 WR stars and LBVs 



iBrevsacher et al.l (|l999l) published a catalogue of 
Population-I Wolf-Rayet (WR) stars in the LMC. 
Among 134 stars, gas mass-loss rates arc given for 40, 
most of which were measured by Crowther fc SmithI 
(|l997t ). Simply integrating these 40 stars gives 2 x 
lO~^M0yr~^. If the remaining 74 stars have simi- 
lar mass-loss rates on average, we expect about 6 x 
10~^ Mq yr~^ from Population-I WR stars in total. 
There are over a hundred suspected WR stars in the 
LMC, so the gas mass- loss rate from WR stars may be 
underestimated by a factor of up to a few. Late-type 
WR stars (WC7-WC9) prod uce dust, but no such stars 
have been found in the LMC (jSmith et al.lll988l : iMoffatI 

[Mil). 

Luminous Blue Variables (LBVs) could also poten- 
tially contribute_to the gas and dust injected into the 
ISM. IZickgraj (I2006D listed 15 B [e] supergiants in the 
LMC and IVan Genderenl (|2001 h listed 21 S Dor vari- 
ables. However, the mass-loss rates of these stars are 
still unknown. Well studied object s of this type in the 
Galaxy include rj Car and AG Car. iHillier et al] (|2001l ) 
suggested a gas mass-loss rat e of 1 x 10~^ Mp, yr~^ fo r 
77 Car using a formula from IWright fc Barlow! (|1975D : 
near-infrared interfe rometric observations suggested 
1.6 X IQ-^MfDyr" ^ (jvan Boekel et al.ll2003f) for 77 Car. 
IVoors et al.l (j2000f) estimated a dust mass-loss rate of 
3 X 10~^ Mq yr~^ for AG Car from mid- and far-infrared 
observations. However, it is clear that not all LBVs have 
dust mass-loss rates as high as 77 Car and fewer than half 
of th e objects have noticeable excesses (IVan Genderenl 
[200lh . 

Winds from OB-type stars co uld contribute to the 
gas enrichment process in the ISM. IBastian et al.l (|2009t) 
estimated about 2000 OB stars to be present within the 
inner 3 degree radius of the LMC. This area is simi- 
lar in size to the SAGE survey coverage. If OB stars 



have mass-loss rate s up to 1 x 10 ^-1 x 10 ^ Mq yr ^ 
(jEvans et al.ll200^ . the total gas lost from OB stars 
could be 2 X 10-'*-2 x lO'^ Mq yr'^. 



4.2.3 Red supergiants 

High mass stars (more than 8Mq) can produce 
dust grains during their red supergiant phase (e.g. 
IVerhoelst et al.|[2009() . One of the most luminous red 
supergiants (RSGs) in the LMC is WOH G64 and its 
circumstellar dust m ass could be as high as 2 x 10^^ Mq 
(jOhnaka et al.ll2008l ). This exceeds the total dust pro- 
duction from SN 1987A. lvan Loon et al.l (|2005l ) listed 14 
oxygen-rich red supergiants, which exceed the classical 
upp er limit of AGB lu minosity, i.e. logL(LQ) = 4.74 
(c.f. IWood et all I1983D . The combined gas mass-loss 
rates of these 14 stars is 1 x 10"'^ Mq yr~ ^ . A gas-to-dust 
ratio of 500 was assumed in their work ()van Loon et al.l 
I1999D . The dust mass-loss rate from RSGs amounts to 
2 X lO-'^Moyr-i. 

Previous studies mainly depend on the IRAS sur- 
vey. After an investigation of RSGs using SAGE data 
(Kemper et al. in preparation), the number of known 
RSGs may increase. The integrated mass-loss rate for 
the RSGs quoted here is therefore likely to be a lower 
limit. 



4.2.4 Oxygen-rich AGB stars 

The total number of oxygen-rich AGB stars with 
circumstellar envelopes is unknown. Thus, their 
total mass-loss rate remains essentially unknown. 
Ivan Loon et al.l (|2005l ) discuss 7 oxygen-rich AGB stars 
and from these we obtain an integrated mass-loss rate 
of 2 X lO"'' Mq yr~i in gas and 0.4 x 10"'^ Mq yr"^ in 
dust, assun iing a gas-to-dust rati o of 500 as in the orig- 
inal work (jvan Loon et al.l |1999| ) . The gas-to-dust ra- 
tio of oxygen-rich AGB stars is unknown, but as silicon 
atoms are not synthesised in AGB stars, the silicate dust 
amount can be assumed to correlate with the original 
silicate abundance of the star. It is expected that the 
gas-to-dust ratios of oxygen-rich AGB stars are lower 
in the LMC than in our Galaxy, as a similar trend was 
found in the ISM. These provide a lower limit to the in- 
tegrated mass-loss rates; further estimates can be made 
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after the IRS foUowup of SAGE is completed (Kemper 
in preparation). 

We expect that the upper hmits of gas and dust 
input from oxygen-rich stars are comparable to the val- 
ues from carbon-rich stars, simply from the relative 
numbers of objects. It is known that the number ratio 
of oxygen-rich st ars and carbon-rich stars depends on 
the metallicities (|Feastlll989l 120061: lGroenewege"^l2007l 
e.g.) In the LMC, among AGB stars with low mass- 
loss ra tes there are as ma ny oxygen-rich as carbon-rich 
stars ([Blanco et al.l llQTSh . However, among PNe, the 
next evolutionary phase, there are on ly two oxygen-rich 
dusty sources found out of 23 PNc (S tanghellini et all 
l2007t iBernard-Salas et all [2008) . Thus, the integrated 
mass-loss rate from oxygen-rich AGB stars is ex- 
pected to be similar to or less than that from carbon- 
rich stars. Further, stellar evolutionary models pre- 
dict that oxygen-rich stars evolve from both low- 
mass (1.0-1.5 M0 on the main sequence) and inter- 
mediate mass (5.0-8.0 Mq) progenitors, while carbon 
stars originate from 1.5-5.0 at the LMC metal- 
licity ([Karakas fc Lattanzioll2007[ ). Many intermediate- 
mass oxygen-ri ch stars will have exp erienced hot bot- 
tom burning ([Boothrovd et al.l I1993D. tu rning carbon 
atoms into oxygen atoms. Kroupa " (|200ll )'s IMF gives 
an oxygen- and carbon-rich ratio of 1:1 in number, 
but 1:1.4 in mass. Hence our conclusion is that the 
integrated mass-loss rate from oxygen-rich AGB stars 
is expected to be similar to or less (approximately 
6 X 10~'^ Mp,) than that from carbon-rich AGB stars. 

iLagadec &: Ziilstral ([2008 al) argue that at low metal- 
licity, oxygen-rich dust drives a superwind only when 
the stellar luminosity exceeds a crictical value, of the 
order of 10^ L© for the case of the LMC. This would 
predict that the mass loss is more dominated by carbon 
stars than their number ratio would suggest. 

Among AGB stars, S-type stars represent stars with 
a C/0 ratio close to unity. This type of star are con- 
sidered to trace a transiting phase from oxygen-rich 
stars to carbon-rich stars. In the LMC, there is no com- 
plete sensu s of S- type stars; a few stars are listed in 
iCioni et all ([200l . In our Galaxy, the fraction of S- 
type stars is assumed to be 1 % of total AGB pop- 
ulation. Thus, we assume gas and dust lost from S- 
type stars is relatively minor. Except for a few SC-type 
stars ([Aoki et al.lll998[ ). which have a carbon- rich atmo- 
sphere, S-type stars tend to p r oduce some silicate d ust 
([Little-Marenin fc Littlel[l988l : iChen fc Kwokl[l99l . 



4.3 Gas and dust sources in the LMC 

Table[ni summarizes the estimated gas and dust produc- 
tion from various sources. AGB stars are one of the main 
sources of the dust enrichment for the ISM of the LMC, 
and carbon-rich AGB stars are a major factor. The 
dust contribution from SNe is very uncertain. There 
could also be minor inputs from novae and RCrB stars. 
Dust evolution models show that within our Galaxy car- 
bon dust originates mainly from carb on-rich AGB stars 
([Dwekl 119981: IZhukovska et al.ll2008l ). Our work shows 
that dust production in the LMC follows a similar trend. 
If a gas-to-dust ratio of 200 is valid for carbon-rich 



AGB stars, and our estimates for mass loss from SNe are 
correct, AGB stars are as important a gas contributor as 
are SNe. SN activity is concentrated towards the 30 Dor 
and LMC bar region. In the LMC disk, AGB stars could 
be the more important source of gas and dust. Other 
sources appear to present minor contributions to the 
ISM enrichment. 

This result is in marked contrast to our Galaxy 
where AGB st ars are thought to be the main source, 
even for gas ([Tielens et al.l l2005l) . The contribution 
from SNe is an order of magnitude smaller than that 
of AGB stars. The observational difference can be re- 
lated t o the recent increase in th e star formation rate 
(SFR; ISmecker-Hane et al.|[2003 ) in the LMC bar re- 
gion, which has been sugge sted to be due to the tidal 
interaction with the SMC ([Bekki fc Chibal l2005f ) . In- 
deed, assuming a constant SFR and using the IMF, we 
find that the gas injection rate from AGB stars should 
be a factor of twice larger than that from SNe and their 
precursors. Thus an increase of the SFR appears to be 
required in the LMC. 

4.4 Populations and spatial distributions 

Table[3| shows, as a function of mass- loss rate, the num- 
ber distributions for carbon-rich AGB candidates in the 
bar, the bar centre, outside of the bar, and the entire 
area surveyed. We have defined the bar centre and bar 
regions as shown in Fig. [51 The fraction of stars with dif- 
ferent mass-loss rates shows little dependence on their 
location within the LMC. This suggests that most of 
these stars belong to a similar population. Contamina- 
tion by non-LMC objects, i.e., distant galaxies, should 
be minimal; the number density of distant galaxies is 
almost uniform and thus the number counted should be 
pro portional to the area examined. 

iReid fc Parked ([2006D surveyed PNe in the LMC us- 
ing narrow-band optical filters. They found fewer PNe 
in the north part of the bar (between 5h45m, — 68°30' 
and 5hl5m, — 67°15') than in the surrounding area, and 
concluded that interstellar extinction towards this re- 
gion is responsible for the lower number density. How- 
ever, a similar trend is found for the AGB-star distri- 
bution, despite the fact that extinction in the near- and 
mid-infrared will be negligible. Generally, in compari- 
son to the age of galaxies, AGB stars and PNe belong 
to similar stellar populations. Thus, this absence of PNe 
and AGB stars in the north part of the bar appears to 
be intrinsic to the LMC. 

A distinct component of the popula tion is also 
found a mong young red supergiants; fig. 2 in lCioni et al.l 
(j2000af) shows a secondary peak at about Dec > —67° 
which presumably corresponds to stars younger than 
0.5 Gyrs. Similarly, high mass-loss rate AGB stars (rates 
between 6 x lO^^Mgyr"! and 3 x 10'^ M^yi-^) 
are populous in that region, although they are ab- 
sent in the northern part of the bar. In this region, 
stars with low and intermediate mass-loss rates (from 
3 X lO-'^Mgyr^i to 6 x lO-'^Mgyr-i) are under- 
populated. This may suggest that the location of active 
star forming regions has not changed significantly over 
the past 0.5-1 Gyrs. 
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Oxygen-rich AGB stars originate both from rel- 
atively low mass (1.0-1.5 Mq) and relatively high 
mas s (5.0-8.0 M(7^) stars, at the LMC metallic- 
ity (jKarakas fc Lattanziol l2007t ). The distribution of 
oxygen-rich and carbon-rich stars might differ from each 
other, as the majority of the oxygen-rich stars (which 
belong to a lowe r mass group) , belon g to a n older pop- 
ulation. Indeed, iBlanco et al.l (|1978D and ICioni et all 
(|2006aD showed that oxygen-rich stars spread further 
out into the disk region than do carbon-rich stars. 

Based on optical imaging an d a population synthe- 
sis model, ISmecker-Hane et al] ([2002) suggested that 
in the bar region 35 percent of the stars are younger 
than 3Gyr, and 71 percent are younger than 7.5 Gyr. 
In contrast, the disk has an older population, as the 
corresponding numbers are 19 percent and 41 percent. 
The initial masses of carbon-rich stars are 1.5- 5 . Mq 
(IVassihadis fc WooJlQQ.lHStancliffe et al.ll2005r) . These 
stars reach the thermal pulsating AGB phase in under 
3 Gyr. Our observations show that the carbon-rich AGB 
candidates are concentrated in the bar, which is consis- 
tent with a population younger than 3 Gyr in the LMC. 



4.5 Evolution of the ISM 

4-5.1 Star formation rate and stellar feedback 

The LMC SFR is higher than the gas injection rate 
from SNe and AGB s tars (in total 0.03-0.07Mo yr'^). 
IWhitnev et all (I2008D estimated a SFR of 0.19 Mq yr"' 



and iKennicutt et al.l (|l995f ) suggested a SFR of 



0.26M©yr~^. These are an order of magnitude higher 
than the gas output from AGB candidates and SNe. 
This suggests that the current LMC star formation de- 
pends on the large rese r voir o f ISM gas (7 x 10^ Mq 
in Hi: IWesterlun d et al.1 (|l997l ) and 1 x 10^ Mq in Ha; 
llsraell (1997)). This is summarised in Table[i With- 
out continuing infall, the LMC would become gas 
poorer, due to the excess SFR, and the SFR in the 
LMC would eventually decline. The Magellanic stream 
has in fact removed further material from the LMC. 
Further, stellar feed back for at least 1 Gyrs is re- 
quired to gain 10% of gas mass present in the ISM, 
suggesting stellar yields can impact the ISM abun- 
dance only over a time scale of a few Gyrs. This 
is consistent with a time scale of the metallicity in- 
creases (order of a Gyr ), accordin g to the age-and- 
meta l licity relation (e.g. Feast 19891: Russell fc Dopital 



119921: iPagel fc Tautvai§ienelll998HCarrera et al.ll2008D . 



4-5.2 Impact on extinction curve 

The main sources of gas and dust are different in our 
Galaxy and in the LMC. This begs two questions: 1) 
why does the gas-to-dust ratio in the ISM in the LMC 
differ from that in the Galaxy? 2) Why do the propor- 
tions of amorphous silicate and carbon dust grains in 
AGB stars appear to differ from those in dust grains in 
the ISM. 

The different enric hment may be expect ed to have 
an impact on the ISM. iGordon et al.l (|2003[ ) measured 
that the gas-to-dust ratio varies within the LMC. 



The highest values was found in the LMC 'super- 
shell' (N(Hl) /A{V)^b-2i X 10^1 atoms cm^^ mag^^) , 
but overall gas-to-dust ratio remains simi- 
lar (N(Hl) /A{V)=l~7xW^ atoms cm'^ mag^^) 
but slightly higher than the Galacitic value 
(N(Hl)M(y)=1.5xl O^^ atoms cm-2 mag-^). 
iGordon et al.l (j2003l ) explained the variation with 
the local condition. The LMC 'supershell' is strongly 
affected by SNe, destroying dust grains. 

The extinction curve in the LMC is quite similar to 
the Galactic one, except for a slight steep rise i n far- 
UV region in 30 Dor region (IGordon et al.ll2003fl . The 
2175 A bump is assumed to be associated with carbona- 
ceous dust grains; this feature is weak in the 30 Dor 
region. Although the LMC has a higher carbonaceous 
dust input for AGB stars, compared to silicate dust, 
than does our Gal axy, this does not ap pear in the ex- 
tinction curve. As IGordon et al.l (|2003D suggested, the 
extinction curve in 30 Dor region is strongly affected by 
SN dust destruction. However, overall the consistency 
of the extinction curve in the LMC and our Galaxy can 
not be explained with SN dust destruction. An addi- 
tional silicate dust source may be required. 



4-5.3 AGB feedback with respect to past SFR in the LMC 
and our Galaxy 

We can estimate the gas and dust feedback rate from 
carbon-rich AGB stars and compare it with the mass 
removed from the ISM by star form ation. We adopt the 
SFR i n the bar and the disk froni iSmecker-Hane et al.1 
(120021) of O.OlM0yr-ideg-2 and 0.004 Mq yr"! deg-^, 
respectively. These numbers approximately trace the 
SFR 1-3 Gyr ago, when most carbon-rich AGB stars 
were formed. The gas ejection rate (GERcagb) and dust 
ejection rate (DERcagb) from carbon-rich AGB stars 
are taken from Table|3| In the bar area, GERcagb/SFR 
is 0.03-0.04 (up to 0.04) in mass, and outside of the bar 
the ratio is 0.02, i.e., the two regions have similar rates. 
This ratio could depend on the IMF, but no apparent 
difference is found in these two GERcagb/SFR. In the 
LMC, more high mass-loss rate AGB stars tend to be 
carbon-rich than oxygen-rich. 

If we include the oxygen-rich stars, GERagb/SFR 
could be higher than GERcagb/SFR by a factor of 1.4- 
2, resulting in GERagb/SFR of 0.03-0.08. If the SFR 
is approximately constant over the life time of the AGB 
stars, about 3-8% of the stellar mass is returned to 
the ISM as gas through AGB stars. The uncertainties 
do not include the effect of the unknown condensation 
temperature of carbon-rich dust grains, which results 
in a systematic error of up to a factor of two on the 
colour and mass-loss rate relation. The ratio of the dust 
ejection rate from carbon-rich AGB stars to the SFR is 
1-2 X 10-*. 

Similarly, we estimate the ratio of the gas 
ejection rate from AGB stars (GERagb) to the 
SFR for Galactic stars, where we obtain 0.02- 
0.06 for GERagb/SFR and 0.8-3x10"! Here we 

adopt the integrated AGB mass-loss ra te of 3- 

exlO^-^Moyr^ikpc-^ (jJura fc KleinmannI [19891 ). In 
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Table 6. Evolution of the ISM 



Gas budget 

Stellar feedback (AGB + SNe) 0.03 - 0.07 Mq yr"! 

Star formation rate (SFR) 0.19 - 0.26 Mq yr"! 

Gas consumption rate (AGB + SNe - SFR) -0.23 0.12 Moyr"! 

Existing gas in the ISM (HI + H2) 8 X 10* Mq 

Time to exhaust ISM gas '-^ 3 X 10^ years 

Dust budget 

AGB dust production rate (5.5-11.5) X 10^^ M0yr~i 

SNo dust production rate (0.01-13) X 10~^ Mq yr'^ 

SNe dust destruction rate ? 

Life time of dust 2-4x10* years 

Accumulated dust mass over the life time of dust (with constant rate) 5 X 10* Mq 

Existing dust in the ISM 1.6 X lO'' Mq 

Age of the LMC 10-15 Gyrs 



this case the sample is not hmited to carbon-rich stars, 
as oxygen-rich stars are included. The SFR of our 
Galaxy is ve ry uncertain and we adopt the value from 
iRanal (|l99lD . who suggest 0.01-0.02 yr"! kpc'^ at 
3 Gyrs ago. The measured AGB feedback rates in the 
LMC and the Galaxy are similar. 

The AGB feedback rate can be compared with the- 
oretical estimates. Stars within the mass range 1.0- 
8.0 Mq comprise abo ut 24 p e rcent of the total mass 
of stars, according to iKroupal (|200lD . If these low- and 
intermediate-mass stars lose about 50-80 percent of 
their mass during the AGB phase, as estimated from the 
white dwarf mass, then 12-19 percent of the feedback 
to the ISM should theoretically be from these stars. 
The AGB feedback rate is only 3-8 percent of the stel- 
lar mass, with a potential increase up to 16 percent. 
In general, our observations are consistent with stellar 
feedback through AGB stars but it appears that the 
theoretical stellar feedback is slightly higher than the 
observed one. 



pected dust mass (5 x IQ'^ Mq) over constant dust sup- 
ply from SNe and AGB stars over 4 x lO^years. This 
shows a 'missing dust-source problem' in the LMC. 

A similar missing dust so urce problem is also 
found in high-z galax ie s (iHughes et al. I 119981 : 
iMorean fc Edmundsl 120031) . (|Meikle et al.ll2007D found 
that the total dust produced from SNe should be far 
less than the dust detected in the high-z g alaxies (e.g. 
IO^Mq; iBertoldi et all I2003D . Recentlv [Sloan et al.l 
(120091 ) found that AGB stars' can contribute to dust 
production in addition to SNe, even for galaxies at 
z~6. But although the dust production AGB stars can 
alleviate the discrepancy, it is unlikely to resolve it in 
high-z galaxies. Additional dust sources are therefore 
needed in the ISM in both the LMC and high-z galaxies. 
On e possibility is star fo r mation regions, as suggested 
by [Morgan fc Edmunds! (|2003D . Alternatively, dust 
processing and dust growth by gas accretion in the ISM 
may yield a total dust mass in the ISM which exceeds 
the dust mass produced by stars. 



4.5.4 Missing dust-source problem 

The total dust mass in the ISM is approximately 1.6 x 
10^ M0. For this estimate, we use the total ISM gas 
mass in the LMC of 8 x 10* Mq, and assume that the 
dust to gas extinction ratio in the LMC i s about a fac- 
tor o f 2-2.5 lower tha n the Milky Way (jGordon et al.l 
l2003t ICox et aP l2006f ) . the adopted gas-to-dust ratio, 
of 500 in the LMC , is higher than the Galactic value 
(jZubko et al.ll2004l ). 

The total dust injection rate from carbon stars is 
measured to be 4.3 x lO~^M0yr~^ and up to 1.0 x 
lO~'*M0yr~^. The rate from oxygen-rich stars is ex- 
pected to be 1.2 X 10~^ Mq yr~^. This is summarised in 
TablelHl Ignoring dust destruction, the dust mass in the 
ISM corresponds to about 20 Gyrs at the current injec- 
tion rate, which is far too long. In fact, the life time of 
dust grains is of the ord er of 2-4 x 10* year s, due to de- 
struction by SN shocks (jJones et al.lll994D . rather less 
than 20 Gyrs. Even adopting the highest dust produc- 
tion rate from SNe, the existing dust mass (8 x 10* Mq) 
in the ISM is an order of magnitude high than the cx- 



5 CONCLUSIONS 

We have analysed Spitzer data for AGB stars in the 
LMC, and presented a quantitative analysis of gas and 
dust budge in the LMC ISM. We obtain the following 
main conclusions: 

• The mass-loss rates of the AGB stars are a quan- 
tifiable function of their infrared colours. 

• The total mass-loss rate from all carbon-rich AGB 
candidates is estimated to be 8.6 x 10^^ Mq yr~^ (up 
to 21 X lO-^Mgyr-i) in gas and 4.3 x lO-^Mgyr"! 
(up to 1 X lO~'*M0yr"i) in dust in the LMC. Adding 
oxygen-rich AGB stars could increase the gas feedback 
rate to 15 x IQ'^ Mq^i'^ (up to 27 x lO-^Mgyr-i) 
and the dust feedback rate to 6 x lO^'"' Mq yr^^ (up to 
1 X lO-^Moyr"^) from AGB stars. 

• AGB stars are among the main sources of gas and 
dust in the LMC. AGB stars are the dominant source of 
gas in the disc, but SNe are a slightly more important 
gas source in the bar region of the LMC. The gas injec- 
tion rate from SNe could be as high as 4 x 10~^ Mq yr~^ 
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in the LMC, but the dust production of SNc is uncer- 
tain. 

• The measured feedback from carbon-rich AGB 
stars into the ISM is 2-4 percent of the SFR 3 Gyr 

ago for gas mass, and 0.01-0.02 percent of the SFR for 
dust mass. The gas feedback from both carbon-rich and 
oxygen-rich AGB stars is expected to be 3 8 percent of 
the SFR. These values could increase by a factor of two, 
due to the uncertainty of dust condensation tempera- 
ture. These ratios appear to be similar in number within 
the LMC bar, the disk and Galactic disk. 

• The gas output of SNe is found to be more impor- 
tant in the LMC bar than in both the remaining region 
of the LMC and in the Galaxy, related to a higher gas- 
to-dust ratio in the LMC. 

• The total gas mass supplied through AGB stars 
and SNe is far less than the current star formation rate. 
The star formation largely depends on the gas already 
present in the LMC ISM, and without infall, eventu- 
ally the SFR will decline as a consequence of the gas 
exhaustion in the ISM. 

• In the LMC, there is 'a missing dust source' prob- 
lem present. The largest estimate of diist accumiilation 
from AGB stars and SNc can not account for the full 
dust mass in the ISM. This is a similar problem as found 
for high-z galaxies. Additional dust sources are required, 
possibly located in star-forming regions. 
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APPENDIX A: CLASSIFICATION OF 
CARBON-RICH STARS 

Al Existing classification schemes 

Fig. lAll shows a J — Ks vs Ks colour magnitude dia- 
gram for the LMC stars from the SAGE catalogue. The 
coding of the chemical types of the stars is given in the 
figure. There are five branches within 11 > Ks > 7. 
Carbon-rich stars are found in the reddest of these 
branche s (J — Ks ^ 1.4), as is kn own from previous 
studies ()Nikolaev fc Weinber gj [2000l ) . From J - Ks ^ 
1.4 to 2.0 mag there is a bend in the carbon-rich branch 
and the ifs-mag becomes brighter with colour. This 
branch reaches its brightest iiTs-mag at about J — Ks ~ 
2, after which Ks becomes fainter. The number den- 
sity also decreases redward of J — Ks ~ 2. Stars with 
J — Ks > 2 are desc r ibed as obscured AGB stars by 
iNikolaev fc Weinberg) (|200ClD . and in our sample most 
of these are carbon-rich, but oxygen-rich st ars are also 
found here. Following the classification by ICioni et al.l 
(|2006a^ {Ks < -0.48 x (J - Ks) + 13 and Ks < 
-13.333 X (J - Ks) + 24.666), only 17 percent of the 
stars are redder than J — Ks > 2. This shows that 
less than one-sixth of carbon-stars have a circumstellar 
envelope detectable in the Ks-ha, nd. This is consisten t 
with the stellar evolution model bv! Marigo et al.l ()2007f ) , 
who showed that the super-wind phase (with mass-loss 
rate higher than 5.6 x 10^^ Mq yr~^) lasts less than 15 
percent of the thermal-pulsing lifetime for carbon-stars. 

The 2MASS sensitivity is not sufficiently high to 
detect dusty carbon-rich AGB stars in some bands. 
Figs.ES and E3] show that 2MASS does not detect the 
J- and ii's-bands at the red end of the carbon-rich se- 
quence. However, this red end is clearly detected at 
[5.8] and [8.0] (Fig.g] and Fig.El). The situation will 
change once the VISTA LMC survey of the Magellanic 
System is cor npleted, which co uld detect as faint as 
ivrs=20.3mag (jCioni et al.ll2008D . 

Oxygen-rich stars are mainly found at Ks^S mag 
and J — Ks < 1.5. These stars a re a mixture 
of re d giants and re d -supe r giants (jKastner et al.l 
|2008[) . ICioni fc Habind (|2003l ) show that oxygen-rich 
AGB stars are found in the branch located at about 
Ks ^9 mag and J — Ks < 1.5. In our analysis fFig. lAl[) . 
there are two regions with oxygen-rich stars . One is the 
same as identified bv lCioni fc Habind (|2003fl . The other 
is occupied by oxygen-rich stars with think dust shell, 
which are just above the obscured carbon-rich sequence. 
There are two oxygen-rich stars with lower luminos- 
ity than the others {Ks ^13mag; J — i4rs=l. 5-2.0): 
163-31 (SAGE 045433.86-692036.1) and MSXLMX610. 
Their low luminosities and some excess at near-infrared 
wavelengths suggest that they could be oxygen-rich 
post-AGB candidates or YSO candidates. Post-AGB 
stars with optica l variability show simi lar J — Ks and 
ifs-magnitudes (|Wood fc CohenI l2001h . The possibil- 
ity that they are YSOs is not discarded, but if so the 
Ks ~ 13 m a.g implies a relatively high mass (between 5 
and 10M,T,: IWhitnev et al.|[200l . 

Fig.|M] shows the J - [3 6] vs [3.6 ] colou r 
magnitude diagram, used by iBlum "etall (l200l . 
in which the separation of oxygen- and carbon- 
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Figure Al. 2MASS J — Ks vs Ks colour magnitude diagram for 
the LMC, extracted from the SAGE catalogue. The lower graph 
shows cross-identifications with carbon-rich AGB stars, oxygen- 
rich AGB stars and RSGs. Large blue circles: carbon-rich AGB 
stars according to the ir IRS spectra. Smal l blue circles: carbon- 
rich AGB stars from jKontizas et al.ll200lh . classified from opti- 
cal spectra. Red triangles: oxygen-rich red giants and supergiants 
based on IRS spectra. 



rich AGB stars is clear. Oxygen-rich stars have 
OH in the J-band, and mos t Mira variables have 
H2O in JHK- and L-bands dlancon fc Woodl I2OOOI: 
i Tsuji et al. 1 Il997l : iMatsuura et all l2002at iTei et al.l 
,2003h . Car bon-rich stars have HCN and C2H2 bands 
at 3.1 um dRidgway et all Il978t iGroenewegenI 11994 



at d.ijum llKidgway et al.l lliJY8t lUroenewegeni iij!j4t 
Ivan Loon et al.lll999HMatsuura et al."2 002bl.l2005D and 



CN in the J- and fl"-bands ([Whitcloc k fc Catchpole I 
119871: iLanpon fc Woodir2000D . The sensitivity is limited 
by the J-band, and heavily dust-obscured AGB stars 
are not detect ed. In Fig.lAll o p tically identified carbon 
stars listed bv iKontizas et all (|2001 ) reach the bright- 
est Ks mag. The Ks magnitudes become fainter as the 
stars become redder in J — Ks, as they do among the 
more obscured carbon stars classified by Spitzer IRS. In 
contrast, such a trend is not found in Fig. lA2[ the [3.6] 
magnitude continues to brighten among obscured car- 
bon stars beyond the J-band detection limit, as would 
be expected. 



A2 Classification using Ks — [8.0] vs Ks 



Fig. lA3l is a colour-magnitude diagram showing Ks — 
[8.0] vs Ks. In this carbon-rich stars become fainter in 
Ks as circumstellar reddening changes the colour; Ks 
reaches a peak of Ks ^ 10 mag around Ks — [8.0] = 2- 
4 mag. Oxygen-rich stars have slightly higher luminosity 
in Ks at a given colour, except for the two oxygen-rich 
post-AGB/YSO candidates. In this diagram, oxygen- 
rich stars are well separated from carbon-rich stars. 
LMC carbon-ri ch stars have HCN and C2H2 absorption 
bands at 7 urn (lAoki et al.lll998l : I jorgensen et al.ll2000l : 
IMatsuura et al.l 120061 ) hTthe [8.0]-band (6.5-9.5 //m) 
while oxygen-rich stars have only weak SiO absorp- 
tion and silicate excess in the [8.0]-band, until the sil- 
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Figure A2. J-band and J — [3.6] colour magnitude diagram. 

icate band turns to absorption. These spectral differ- 
ences separate tlie two chemical types in this colour- 
magnitude diagram. This diagram also shows that high 
mass-loss rate carbon-rich stars could have Ks-h&nd 
magnitudes fainter than 15 mag, i.e., below the detec- 
tion limit of the 2MASS A's-band. 

A3 Other colour-magnitude diagrams 

High mass-loss rate stars are likely to be detected at 
longer wavelength, as is illustrated in the Spitzer SAGE 
colour-magnitude diagrams discussed below. Figs. lA4l 
and IA5I show colour-magnitude diagrams using the 
SAGE [5.8]-, [8.0]- and [24]-bands. In these two figures 
the sequence of IRS-classified carbon- rich stars does not 
end at the detection limits. Thus all of the red, i.e., high 
mass-loss rate carbon-rich, stars were detected in these 
bands. However, there is no clear separation between 
oxygen-rich and carbon-rich stars so the diagrams are 
not useful for classification. 

A4 Near-infrared survey with the IRSF 

iKato et ai] (|2007D pubhshed a JHKs survey of the 
Magellanic Clouds taken with the Infrared Survey Facil- 
ity (IRSF) telescope at SAAO in South Africa. The sen- 
sitivity of this survey, 16.6 mag at 10 cr in the ifs band 
for the LMC, is significantly better than that of 2MASS 
(14.3 mag), although the conversion between IRSF and 
2MASS filt er system may n ot be straightforward for 



A GB stars (iKato et al. | l2007D . Four of the stars hsted 
in iGroenewegen et al.l ( 2007^ have been observed with 
IRSF and their mass-loss rate is plotted in Fig. lA6l It 



appears that IRSF stars follow a similar colour mass- 
loss rate trend. We will investigate this relation further 
after the Spitzer spectroscopic survey, SAGE-Spec, is 
complete. 
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Figure A3. Ks and Ks — [8.0] colour magnitude diagram. The 
upper panel shows all of the SAGE sample, and spectroscopically 
known oxygen-rich AGB/RSGs, while carbon-rich AGB stars are 
indicated in the lower panel. The /fs-band sets the detection limit 
at Ks — [8.0] < 2 mag, except for blue {Ks — [8.0] < 2 mag) stars, 
which are [8.0]-band limited. The number of carbon-rich AGB 
stars decreases at about Ks — [8.0]=1.4 
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Figure A4. [5.8] — [8.0] vs [8.0] colour-magnitude diagram. Sym- 
bols are the same as in Fig lAll Oxygen-rich and carbon-rich stars 
are found close together. 



03 

Cl 

cn 
M 

E 



(0 



1 

0} 



10° 



10 



ID' 




2 4 6 

Ks-[8.0] 

Figure A6. The iCs-[8.0] using IRSF Ks is plotted 
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Figure A5. [8.0] — [24] vs [8.0] colour-magnitude diagram. Sym- 
bols are the same as Fig lAll Oxygen-rich and carbon-rich stars 
are not separated. 



Figure A7. [3.6] — [24] vs [24] colour-magnitude diagram. Sym- 
bols are the same as lAll 



